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Abstract 


Cell-to-cell heterogeneity within an isogenic population has been observed in prokaryotic 
and eukaryotic cells. Such heterogeneity often manifests at the level of individual protein 
abundance and may have evolutionary benefits, especially for organisms in fluctuating envi- 
ronments. Although general features and the origins of cellular noise have been revealed, 
details of the molecular pathways underlying noise regulation remain elusive. Here, we 
used experimental evolution of Saccharomyces cerevisiae to select for mutations that 
increase reporter protein noise. By combining bulk segregant analysis and CRISPR/Cas9- 
based reconstitution, we identified the methyltransferase Hmt1 as a general regulator of 
noise buffering. Hmt1 methylation activity is critical for the evolved phenotype, and we also 
show that two of the Hmt1 methylation targets can suppress noise. Hmt1 functions as an 
environmental sensor to adjust noise levels in response to environmental cues. Moreover, 
Hmt1-mediated noise buffering is conserved in an evolutionarily distant yeast species, sug- 
gesting broad significance of noise regulation. 





Introduction 


Genetically identical cells grown under homogeneous conditions can still exhibit heteroge- 
neous phenotypes. This heterogeneity is ubiquitous and manifests at different levels, from 
individual protein concentrations (protein noise) [1] to cell physiology (cellular noise) [2,3]. 
Although phenotypic heterogeneity only exists transiently, it can lead to deterministic out- 
comes. In multicellular organisms, a stochastic difference in the initial cell state can result in 
different cell fates during development [4,5]. Moreover, stochastic variation in gene expression 
has been shown to determine the outcome of inherited detrimental mutations [6,7], represent- 
ing a possible cause for the incomplete penetrance observed in many human diseases. In 
microbial cells, levels of preexisting heterogeneity can influence population fitness upon expo- 
sure to unpredictable environmental change [8,9]. This “bet-hedging strategy” is commonly 
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used by microorganisms to ensure population survival without the fitness cost of developing 
complex regulatory networks that respond to randomly fluctuating environments [10]. 

At the gene expression level, preexisting cell-to-cell heterogeneity mainly originates from 
the stochasticity inherent to molecular processes (such as transcription factor binding to target 
sequences) and fluctuating levels or activities of factors critical to those processes (such as 
RNA polymerase II or ribosomes) [1,11,12]. Genome-wide studies have shown that low-abun- 
dance proteins often present higher protein noise, which is consistent with the greater variabil- 
ity of infrequent events [13,14]. 

However, some studies have revealed certain pathway-specific patterns in noise levels. For 
example, housekeeping genes tend to have lower protein noise, whereas environment-respon- 
sive genes are often noisier [15,16], perhaps because fluctuations in housekeeping genes may 
compromise essential cellular functions and noisy environment-responsive genes can exert a 
bet-hedging function. The observed patterns in these two types of genes indicate that selection 
operates on protein noise levels or that levels have been adjusted according to potential costs 
and benefits over the course of evolution [17]. Moreover, a study comparing young and old 
mice showed that heart cells isolated from old mice exhibit higher cellular noise than those iso- 
lated from young mice [18], suggesting that noise levels are tightly controlled in young healthy 
cells but that the control systems deteriorate with age. 

How do cells adjust protein noise? Several general features have been associated with pro- 
tein noise, including network topology, cellular compartmentalization, molecular chaperone 
abundance, nucleosome occupancy, and promoter architecture [7,19-22]. Genetic studies 
have also identified mutations that alter local or general noise levels [5,12,23,24]. Nonetheless, 
how cells respond to growth conditions and integrate different pathways to fine-tune protein 
noise remains insufficiently characterized. 

To understand how protein noise is regulated, we used experimental evolution of S. cerevi- 
siae to search for mutations that increased the protein noise of different reporter genes. After 
35 cycles of selection, two of the evolved lines (TDH2-GFP- and TYS1-GFP-carrying lines) 
exhibited increased noise levels without a concomitant reduction in protein abundance. We 
show that increased protein noise in the evolved line carrying the TDH2 reporter gene is not 
specific to the TDH2-related pathway, suggesting that the evolved mutations have a general 
effect on protein noise regulation. We identified the methyltransferase Hmt]1 as the major con- 
tributor of the evolved phenotype. Further experiments revealed that noise regulation is medi- 
ated by methylation of multiple downstream targets of Hmt1 and that HMT1 expression is 
often attenuated under stress conditions. Our results suggest that Hmt1 functions as a master 
regulator that adjusts noise levels in response to environmental stress. 


Results 
Experimental evolution of increased protein noise in budding yeast 


A previous study showed that alternating selection between highest- and lowest-expression 
subpopulations could efficiently enrich promoter variants for high transcriptional noise in 
bacterial cells [25]. We hypothesized that a similar selection strategy might allow us to “evolve” 
yeast cells to increase the protein noise of reporter genes (Fig 1A). Eight genes (ADK1, APA1, 
PCM1, RPL4B, SAM4, TDH2, TPD3, and TYS1) selected from distinct cellular pathways were 
fused with green fluorescent protein (GFP) to generate our reporters (see Materials and meth- 
ods). Evolving lines carrying individual reporter genes were subjected to alternating selection 
between the top 5% and bottom 5% of total populations in terms of their GFP intensity. We 
also treated cells with a mutagen (2.8% ethyl methanesulfonate [EMS]) before each selection 
cycle to increase the genetic diversity of evolving populations. After 35 cycles of selection, half 
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Fig 1. Yeast cells exhibit high noise levels upon experimental evolution. (A) Schematic of our evolution experiment. Cells were 
treated with 2.8% EMS, regrown for 12 h and then selected for the top 5% (or bottom 5%) of the population in terms of GFP intensity. 
Alternating selection enriched for mutant cells that presented higher expression noise in terms of GFP intensity. (B) Most of the 
evolved clones exhibit increased expression noise. Single clones were isolated from eight evolved cultures, and their reporter gene 
expression was measured. The x- and y-axes represent fold-change in noise and mean expression, respectively, after evolving. The 
median and range for 4-5 replicates for each evolved clone are indicated by circles and error bars, respectively. The red lines indicate 
values for the ancestral line. Data associated with this figure can be found in the supplemental data file (S1 Data). EMS, ethyl 
methanesulfonate; GFP, green fluorescent protein. 


https://doi.org/10.1371/journal.pbio.3000433.g001 
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of the evolved lines (including the strains carrying RPL4B, TDH2, TPD3, or TYS1 reporter 
genes) exhibited significantly increased reporter noise (Fig 1B). 

Among them, the Tdh2-GFP- and Tys1-GFP-carrying lines also presented increased mean 
protein intensity, ruling out the possibility that the increased noise was due to reduced mean 
protein signal intensity. These results indicate that our selection regime effectively increased 
protein noise upon experimental evolution. In order to dissect the genetic basis of noise regu- 
lation, we selected the Tdh2-GFP-carrying line for further analysis because it exhibited the 
greatest increase in noise without a concomitant decrease in mean signal intensity. 


Isogenic cells with high and low Tdh2 levels have fitness advantages under 
different conditions 


Bet-hedging is a commonly adopted survival strategy among microorganisms for spreading the 
risk of encountering hostile environments [9,26]. Tdh2 protein is a glyceraldehyde-3-phos- 
phate dehydrogenase involved in glycolysis/gluconeogenesis, and it has been shown to help 
cells resist oxidative stress during the stationary phase [27]. We tested whether cells with high 
or low Tdh2-GFP levels represent different physiological states and whether they exhibit fitness 
advantages under different conditions. To do this, we isolated individual stationary-phase cells 
presenting different levels of Tdh2-GFP using a cell sorter and examined their phenotypes. 
Consistent with a previous observation [27], upon HO, treatment, cells with high 
Tdh2-GFP levels had higher survival rates than those with low Tdh2-GFP (Fig 2A). Interest- 
ingly, when we provided fresh nutrients, the cells with high levels of Tdh2-GFP tended to reen- 
ter the cell cycle more quickly than those with low Tdh2-GFP levels, despite cells with different 
levels of Tdh2-GFP exhibiting no difference in survival rates under this condition (Fig 2B). 
This variation in cell cycle reentry is reminiscent of the divergent germination times among 
the individuals derived from the same population in some plants or fungi, which has been sug- 
gested to be a risk-spreading strategy to enhance long-term survival [28-30]. Divergent cell 
cycle reentry times can prevent an entire cell population from going extinct upon occurrence 
of an unpredicted environmental catastrophe. To test this hypothesis, we collected cells dis- 
playing either high or low Tdh2-GFP and challenged the cells with heat stress either before or 
after the cells had been re-fed with fresh nutrients. We found that cells with low Tdh2-GFP 
had a survival rate 3-fold greater than that of cells with high Tdh2-GFP upon encountering 
heat stress after nutrient refreshment. Survival rates were similar and independent of 
Tdh2-GFP levels for cells either nonstressed or stressed before nutrient refreshment (Fig 2C). 
We also found that high and low levels of Tdh2-GFP likely represent transient states, but 
not genetic modifications, of the cells. When we propagated cells sorted into high- and low- 
Tdh2-GFP populations, Tdh2-GFP intensities of both populations reverted to a level similar to 
that of the initial unsorted population after a few generations (S1 Fig). Together, our results 
demonstrate that having cells with different levels of Tdh2 in an isogenic population is advan- 
tageous in different environments, revealing the risk-spreading benefit of expression noise. 


Increased noise is not limited to the TDH2-related pathway in the 
Tdh2-GFP evolved line 


Before characterizing the detailed phenotypes of the evolved Tdh2-GFP-carrying line, we 
examined the cell populations and found that the signal for increased noise was not bimodal 
(S2A Fig). We further confirmed that Tdh2-GFP retained its full length and subcellular locali- 
zation after evolution (S2B and S2C Fig) and that there were no mutations in its promoter or 
coding regions. These data indicate that the increased noise in the evolved line is not caused by 
mutations in the TDH2 locus. 
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Fig 2. Cells with different levels of Tdh2-GFP exhibit different physiological states. (A) Stationary-phase cells with high Tdh2-GFP levels 
(red circles) survive better than those with low levels (blue squares) after growing in H,O -containing medium (Fisher’s exact test, p = 0.1099 
for untreated samples, p < 2.2 x 107'* for H,O,-treated samples). Single cells with high or low GFP intensities were sorted and plated on the 
same plates with or without 4.4 mM H,Ob. Survival rates were determined by counting colony-forming units after 5 d of growth. (B) 
Stationary-phase cells with high Tdh2-GFP tend to reenter the cell cycle faster than those with low Tdh2-GFP signal. Each dot of the 
scatterplot represents data from a single cell. Unsorted stationary-phase cells were placed on YPD agarose pads, and only unbudded cells 
were monitored using time-lapse microscopy. The x-axis indicates the initial Tdh2-GFP signal intensity for each cell, and the y-axis indicates 
the first rebudding time. The red dotted line represents a linear regression (Spearman’s rank correlation, n = 111, p = 8.131 x 10°17). (C) 
Stationary-phase cells with low Tdh2-GFP signal survive better than those with high signal when cells encounter heat stress after being re-fed 
with fresh nutrients (one-sided Wilcoxon rank-sum test, n = 4-5; p = 0.1995 for untreated samples, p = 0.2426 for heat-treated samples before 
nutrient refreshment, p = 0.0159 for heat-treated samples after nutrient refreshment). Survival rates were determined by counting colony- 
forming units after 3 d of growth. The median value for replicates is indicated by solid horizontal lines among groups of data points. 

*p < 0.05, ***p < 0.001. Data associated with this figure can be found in the supplemental data file (S1 Data). GFP, green fluorescent protein; 
YPD, 1% yeast extract, 2% peptone, and 2% dextrose. 


https://doi.org/10.1371/journal.pbio.3000433.g002 


If the evolved mutations in the Tdh2-GFP-carrying line occurred within sequences pertain- 
ing to general noise regulators, these mutations may affect both the TDH2-related pathway 
and other pathways. We selected four pertinent genes—TDH3, PGK1, ADK1, and GLY1—to 
investigate the effect of the evolved mutations. TDH3 is a paralog of TDH2 that encodes an 
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enzyme for glycolysis/gluconeogenesis. PGK1 encodes another key enzyme involved in glycol- 
ysis/gluconeogenesis. Both TDH3 and PGK1 are likely coregulated with TDH2 because all 
three respective proteins operate in the same metabolic pathway [31]. ADK1 encodes an ade- 
nylate kinase required for purine metabolism [32], and GLY1 encodes a threonine aldolase 
involved in glycine biosynthesis [33], neither of which is related to the TDH2 pathway. None- 
theless, all four genes exhibited increased protein noise in the evolved line (Fig 3A), suggesting 
that the effect of the evolved mutations is not pathway-specific. 

We also subjected a TDH2 promoter-driven GFP construct to the same analysis and 
observed increased noise in the respective evolved line (Fig 3B), raising the possibility that 
increased noise in the evolved cell lines may be attributable, at least in part, to promoter regu- 
lation or mRNA turnover. 


An Hmt1 methyltransferase mutant significantly contributes to the evolved 
increase in noise 


To understand the genetic basis of the increased noise we observed in the evolved Tdh2-GFP- 
carrying line, both ancestral and evolved lines were subjected to whole-genome sequencing. A 
total of 1,022 mutations (including 494 nonsynonymous, 256 synonymous, and 271 intergenic 
mutations) were identified in the evolved genome (S1 Table). This high number of mutations 
is most likely due to the mutagen treatment we applied during the cycles of selection. 

Next, we used bulk segregant analysis to refine the list of candidate mutations. To do this, 
we crossed the evolved line to the ancestral line and analyzed their F1 haploid progeny. We 
measured the noise level of Tdh2-GFP in 360 segregants and established an “evolved-like” pool 
(comprising 16 segregants) and an “ancestral-like” pool (20 segregants; see S3 Fig and Materials 
and methods). Both of these segregant pools were then subjected to whole-genome sequencing. 
Based on our computational simulation, we applied two criteria to select candidate mutations 
from these segregant pools: (1) the mutation frequency in the “evolved-like” pool had to be 
>70%; and (2) the difference in mutation frequencies between the “evolved-like” and “ances- 
tral-like” pools was >38% (see Materials and methods). Twenty nonsynonymous mutations 
met these two assumptions and were subjected to reconstitution experiments (S2 Table). 

We first introduced the candidate mutations into the ancestral line using the CRISPR/Cas9 
system and then examined the noise level of Tdh2-GFP. Among the tested reconstitution lines, 
we found that a mutation (G70D) in the methyltransferase-encoding gene, HMT1, resulted in 
a significant increase of noise, close to the level exhibited by the evolved line (Fig 4A). Simi- 
larly, we observed increased noise in the TDH2 promoter-driven GFP construct (Fig 4B). 
When we rescued the mutation in the evolved line by reverting to the wild-type sequence, 
Tdh2-GFP noise was reduced to a level similar to the ancestral line (Fig 4A). These data indi- 
cate that hmt1-G70D is a primary mutation contributing to the evolved increase in noise. 

Already in this study, we have shown that proteins in both TDH2-related and TDH2-unre- 
lated pathways exhibited increased noise in the evolved line (Fig 3A). In the hmt1-G70D recon- 
stitution line, we observed similar noise increases in the TDH3, ADK1, and GLY] reporter 
genes (Fig 4C), indicating that the increased noise observed in the evolved line is mainly due 
to the hmt1-G70D mutation. 


The noise buffering effect of Hmt1 is mediated through multiple 
methylation targets 
Hmtl is a methyltransferase that methylates arginine residues in its substrates. The mutated 


G70D residue is located within a conserved methyltransferase motif (S4A Fig), and mutations 
in this motif cause loss of methyltransferase activity in Escherichia coli and yeast [34,35]. The 
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Fig 3. Multiple genes in the evolved cells exhibit increased noise levels. (A) Four reporter genes (TDH3-BFP, PGK1-BFP, 
ADK1-GEP, and GLY1-BFP) were engineered in modified ancestral (red circles) and evolved (blue squares) TDH2-GFP-carrying cells 
(see Materials and methods), and their expression noise was measured (one-sided Wilcoxon rank-sum test, n = 7-8; p = 4.534 x 10+ 
for ADK1, p = 0.0027 for GLY1, p = 9.441 x 10~* for PGK1, p = 7.158 x 10~* for TDH3). TDH3 and PGK] are involved in TDH2-related 
pathways, whereas ADK1 and GLY] are not. (B) Transcriptional regulation is responsible for increased noise in the evolved cells (one- 
sided Wilcoxon rank-sum test, n = 7-8; p = 4.495 x 10). The TDH2 promoter was directly fused with the coding sequence of GFP. 
This construct was engineered in modified ancestral and evolved TDH2-GFP-carrying cells, and its expression noise was measured. The 
median value of replicates is indicated by horizontal solid lines among groups of data points. **p < 0.01, ***p < 0.001. Data associated 
with this figure can be found in the supplemental data file (S1 Data). BFP, blue fluorescent protein; GFP, green fluorescent protein. 


https://doi.org/10.1371/journal.pbio.3000433.g003 
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Fig 4. The hmt1-G70D mutant recapitulates the increased noise phenotype observed in the evolved line. (A) Reconstituting the 
hmt1-G70D mutation in the ancestral background increases Tdh2-GFP noise, whereas reversing that mutation in the evolved 
background decreases it (one-sided Wilcoxon rank-sum test, = 6; p = 0.0011 for the ancestral background, p = 0.0025 for the 
evolved background). (B) Reconstituted hmt1-G70D mutant cells exhibit increased expression noise of the TDH2 promoter-GFP 
construct (one-sided Wilcoxon rank-sum test, n = 8; p = 0.0002). (C) Reconstituted hmt1-G70D mutant cells present increased noisy 
expression of Adk1-GFP, Gly1-BFP, and Tdh3-BFP (one-sided Wilcoxon rank-sum test, n = 7-8; p = 0.0030 for ADK1, p = 0.0011 for 
GLY1, p = 0.0030 for TDH3). All mutants were constructed using the CRISPR/Cas9 system. The median value of replicates is 
represented by horizontal solid lines among groups of data points. **p < 0.01, ***p < 0.001. Data associated with this figure can be 
found in the supplemental data file (S1 Data). BFP, blue fluorescent protein; GFP, green fluorescent protein. 


https://doi.org/10.1371/journal.pbio.3000433.g004 
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hmt1-G70D mutant also exhibited a noise increase similar to levels exhibited by HMT1-dele- 
tion cells (S4B Fig). Moreover, western blot analysis using anti-methylated arginine antibodies 
confirmed that hmt1-G70D mutant cells lack Hmt1 methyltransferase activity (S4C Fig). These 
results suggest that the noise-regulating function of Hmtl is mediated through its methyltrans- 
ferase activity. 

In vivo and in vitro studies have identified several Hmt1 substrates [36-38]. Moreover, 
Hmtl-mediated methylation has been shown to enhance the function of its substrates in mul- 
tiple pathways, including chromatin remodeling and transcription [39-41], translation and 
ribosome biogenesis [42-45], and posttranscriptional regulation [46-49]. We selected the rep- 
resentative proteins Npl3, Rps2, Sbp1, and Snf2 from among these pathways and generated 
corresponding deletion or hypomorphic mutants (in cases when mutant haploids died or had 
severe growth defects) and investigated noise levels among these mutants. 

Of all these tested mutants, only the rps2 and snf2 mutants presented significantly increased 
noise (Fig 5A), with the other mutants showing either no change or reduced noise (Fig 5B). 
Rps2 is a component of the small ribosomal subunit, suggesting involvement of translational 
regulation in noise control. Snf2 is the catalytic subunit of the SWI/SNF chromatin remodeler, 
and its function depends on two other SWI/SNF components, Snf5 and Snf6 [12,50,51]. We 
further assessed snf5 and snf6 mutants and found that they also presented increased noise (Fig 
5A and 5C), confirming the general role played by the SWI/SNF chromatin remodeler in noise 
regulation. 

To test whether hmt1-G70D mutation affects the activity of SWI/SNF complexes, we per- 
formed a functional assay of the SWI/SNF chromatin remodeler. It has been shown previously 
that transcriptional activation of the CHA] gene is attenuated when cells have defective SWI/ 
SNF complexes [52]. We monitored mRNA levels of CHA1 immediately after shifting cells to a 
CHA1-inducing medium. Expression of CHA1 was indeed reduced in hmt1-G70D mutant 
cells and in snf5 mutants (Fig 5D), indicating that the function of SWI/SNF complexes is com- 
promised in hmt1-G70D mutant cells. 

We then examined the promoter regions of TDH2 and GLY1 (unrelated to TDH2) by 
histone chromatin immunoprecipitation (ChIP) combined with quantitative PCR (Q- 

PCR) to understand how Hmtl influences gene expression. We observed increased nucleo- 
some occupancy in hmt1-G70D mutant cells (Fig 5E), providing evidence for the supposi- 
tion that noise regulation partially operates via SWI/SNF-mediated nucleosome remodeling 
[21]. 


Hmtl1 functions as a mediator in responses to environmental stresses 


Genome-wide studies have indicated that HMT1 expression frequently fluctuates under differ- 
ent growth conditions [53,54], suggesting that cells may employ an environmental sensor to 
enhance phenotypic heterogeneity upon encountering stress. We measured the mRNA levels 
of HMT1 under various stress conditions—including heat, oxidative stress, high osmolarity, 
and glucose starvation—and found that they were significantly reduced under all of these con- 
ditions (Fig 6A). Consistent with the concept of an environmental sensor, the noise level of 
Tdh2-GFP was also increased under stress conditions (Fig 6B and S5 Fig). Together, these data 
indicate that Hmtl can serve as a mediator to control levels of phenotypic heterogeneity in 
response to environmental stimuli. 

Our data show that bet-hedging within a population could be represented by the noisy 
expression of Tdh2-GFP (Fig 2). If increased noise can help populations survive stressful envi- 
ronments, we anticipated that cells with low Hmtl activity under stress would present an 
enhanced survival rate relative to those with normal Hmtl activity under stress. Indeed, when 
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Fig 5. Mutations of SWI/SNF components and the small ribosomal subunit Rps2 result in elevated noise. (A) Attenuating 
RPS2, SNF2, or SNF5 gene expression in the ancestral background results in increased noise (one-sided Wilcoxon rank-sum test, 
n= 6-19; p = 0.0039 for rps2-DAmP/+, p = 0.0025 for snf2A/+, p = 0.0008 for snf5A/+). Rps2 and Snf2 are methylation substrates 
of Hmtl. Snf2, Snf5, and Snf6 are essential components of the SWI/SNF chromatin remodeler. Since haploid mutants of rps2- 
DAmP, snf2A, and snf5A exhibit severe growth defects and because slow growth has been shown to increase expression noise [91], 
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we constructed heterozygous mutant diploids to measure noise. (B) Deletions of NPL3 or SBP1 in haploid cells do not result in 
increased noise (one-sided Wilcoxon rank-sum test, n = 5; p = 0.0079 for npl3A, p = 0.206 for sbp1A). Npl3 and Sbp1 are 
methylation substrates of Hmt1. (C) Deleting SNF6 in the ancestral background results in increased noise (one-sided Wilcoxon 
rank-sum test, n = 5; p = 0.0087). Haploid snf6A mutants exhibit defective respiration, so we used a haploid rho- ancestral strain in 
this experiment, and cells were grown in YPD medium. The median value of replicates is represented by horizontal solid lines 
among groups of data points. (D) The activity of SWI/SNF chromatin-remodeling complexes is compromised in the hmt1-G70D 
mutant. It has been shown previously that compromising SWI/SNF complexes results in attenuated transcriptional activation of 
CHAI [52]. Total RNA was isolated from ancestral, hmt1-G70D, and snf5A haploid cells 0, 5, and 15 min after adding 0.1% 
L-Serine (an inducer of CHA1 expression). Specific mRNA levels were assessed by Q-PCR. In the figure, CHA1 mRNA levels were 
normalized to those of PYK1, with this latter acting as an internal control. For all time points, CHA1 mRNA levels in ancestral cells 
(blue) were significantly higher than those in hmt1-G70D (red) or snf5A (green) mutant cells (one-sided Student's f test, n = 9; 

p < 0.01). (E) hmt1-G70D mutant cells display higher nucleosome occupancy in the promoter and partial coding regions of TDH2 
and GLY1. The chromatin status of the TDH2 promoter was established by ChIP against H3 coupled with Q-PCR. The numbers in 
the x-axis indicate the distance (in bp) from the transcription start site. The y-axis represents relative enrichment of H3 signals for 
the amplicons at the indicated regions (one-sided Student's f test, n = 10-12). Black, gray, and dashed circles indicate confirmed, 
fuzzy, and condition-specific nucleosome-occupied regions, respectively. Error bars represent standards errors. *p < 0.05, 

**p < 0.01, ***p < 0.001. Data associated with this figure can be found in the supplemental data file (S1 Data). ChIP, chromatin 
immunoprecipitation; DAmP, decreased abundance by mRNA perturbation; GFP, green fluorescent protein; H3, histone 3; 
Q-PCR, quantitative PCR; YPD, 1% yeast extract, 2% peptone, and 2% dextrose. 


https://doi.org/10.1371/journal.pbio.3000433.g005 








we treated wild-type and hmt1-G70D mutant cells with HO , the mutant population exhibited 
better viability (Fig 6C). 


Hmt1-mediated noise buffering is conserved in S. pombe 


Our results suggest that the methyltransferase Hmtl may function as a general noise regulator 
that constrains physiological noise in normal conditions but facilitates heterogeneity under 
stress. This buffering mechanism is likely to increase long-term population survival. Methyl- 
transferase is a conserved enzyme that exists in all eukaryotic kingdoms. We examined 
whether its buffering function is also conserved in a phylogenetically distinct species, 

S. pombe, which diverged from the common ancestor of S. cerevisiae at least 400 million years 
ago [55]. We generated yellow fluorescent protein (YFP) fusion protein constructs of tdh1 and 
gpd3 (the S. pombe orthologs of TDH2) and examined their noise levels in the wild-type and 
rmtl (the S. pombe ortholog of HMT1) deletion backgrounds. Similar to our findings for S. cer- 
evisiae evolved lines, noise levels of Tdh1-YFP and Gpd3-YFP were significantly increased in 
rmt1 mutant cells (Fig 6D). Accordingly, the Hmt1-mediated noise buffering system may rep- 
resent an important survival strategy that is conserved across diverse microorganisms. 


Discussion 


Microorganisms frequently face changing environments. Despite being equipped with com- 
plex stress-adaption systems, unpredicted acute stress remains a challenge for cells. Cell popu- 
lations harboring heterogeneous physiological states enhance their likelihood of surviving 
environmental fluctuations [56,57]. Many examples of bet-hedging have been reported previ- 
ously, indicating that this isa common survival strategy among microbes [5, 58]. However, it 
is not known whether cells exhibit another layer of regulation that allows them to adjust their 
levels of noise. Our current study provides evidence that Hmt1 can function as a core regulator 
to constrain or facilitate phenotypic heterogeneity in response to environmental stimuli. 

The evolutionary advantage of an environment-sensing noise regulator is readily conceiv- 
able. Although stochastic noise is inevitable among individuals within a population, excessive 
deviation from “normal” levels may result in a significant reduction of fitness under normal 
conditions. By regulating multiple pathways, Hmt1 allows cells to constrain noise levels. How- 
ever, when a population is exposed to mild stresses, Hmt1 expression is immediately down- 
regulated, and the noise buffering system is curtailed. The resulting enhanced heterogeneity 
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Fig 6. Hmt1-mediated noise suppression can be released under nonoptimal growth conditions as a conserved cell survival strategy. (A) HMT1 
transcripts are down-regulated under nonoptimal growth conditions (one-sided paired Student’s ¢ test, n = 3-4; p = 0.0008 for heat stress, p = 0.003 for 
oxidative stress, p = 0.002 for osmotic stress, p = 0.007 for diauxic shift). The level of mRNA was measured using Q-PCR, and DET1 mRNA was used as 
the internal control. Error bars indicate standard errors. (B) Nonoptimal growth conditions result in increased Tdh2-GFP noise (one-sided Wilcoxon 
rank-sum test, n = 4; p = 0.015 for heat and oxidative and osmotic stress, p = 0.014 for diauxic shift). Noise was measured 2.5 h after cells had been 
shifted to the indicated conditions or was measured at different growth stages. (C) hmt1-G70D mutant populations survive better than the ancestral line 
in medium containing HO? (Fisher’s exact test; p < 2.2 x 101°). Stationary-phase cells were spread on plates with 5.3 mM H20Oz, and survival rates 
were determined by counting colony-forming units after 5 d of growth. (D) Mutation of the HMT1 ortholog in the fission yeast Schizosaccharomyces 
pombe also results in increased expression noise. rmt1 is the ortholog of HMT1, whereas tdh1 and gpd3 are TDH2 orthologs. Deletion of rmt1 increased 
protein noise of both tdh1-YFP and gpd3-YFP (one-sided Wilcoxon rank-sum test, n = 10-12; p = 0.0005 for tdh1-YFP, p = 0.037 for gpd3-YFP). The 
median value of replicates is indicated with horizontal solid lines among groups of data points. *p < 0.05, **p < 0.01, ***p < 0.001. Data associated with 
this figure can be found in the supplemental data file (S1 Data). GFP, green fluorescent protein; Q-PCR, quantitative PCR; YFP, yellow fluorescent 
protein. 


https://doi.org/10.1371/journal.pbio.3000433.g006 


may allow the population to increase its likelihood of survival, especially if the stress is pro- 
longed or escalates (Fig 7). 

Hmtl is a type I arginine methyltransferase that catalyzes monomethylation and asymmet- 
ric dimethylation in budding yeast. Hmt1-mediated methylation has been shown to influence 
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Fig 7. A model showing how Hmt1 modulates cell-to-cell heterogeneity in response to environmental stress. Hmt1 
methylates and enhances the function of the SWI/SNF chromatin remodeler and small ribosomal subunits to reduce 
stochastic noise in gene expression. Under normal conditions, cells maintain a high level of Hmt1 and exhibit 
homogeneous gene expression in most cells. However, when the population encounters environmental stress, HMT1 
expression is down-regulated, inhibiting the functions of Hmtl targets. Accordingly, expression of Hmtl gene targets 
becomes noisier, so individual cells exhibit heterogeneous cell physiologies. The likelihood of population survival is 
enhanced because of this heterogeneity. Me, methylation. 





https://doi.org/10.1371/journal.pbio.3000433.g007 


various cellular pathways, including transcription, posttranscription, and translation [39- 
41,43,45-49]. Moreover, Hmtl is a highly interactive protein, ranking among the top 1.5% of 
the entire budding yeast proteome in terms of the interaction number [59], suggesting it has a 
central role in coordinating a complex network. Previous studies have shown that another net- 
work hub, heat shock protein 90 (Hsp90), can regulate general cellular noise [20,60]. Consis- 
tent with predictions from network analyses and results from a genome-wide screen, hub 
genes have a strong impact on buffering nongenetic and genetic variation [61-63]. More inter- 
estingly, both Hsp90 and Hmt1 can work as environmental sensors to adjust noise levels in 
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response to environmental cues, acting as direct links between noise regulation and adaptive 
benefits. 

Our results show that at least two downstream pathways, i.e., chromatin remodeling and 
the translational machinery, are involved in Hmt1l-mediated noise buffering. The effect of the 
SWI/SNF chromatin-remodeling complex on gene noise has been reported previously [12]. 
Here, we have identified Hmtl as an upstream coordinator of noise regulation and establish 
the novel role of protein methylation in this process. The detailed mechanism underlying how 
the ribosomal component Rps2 regulates noise is less clear. Rps2 plays a crucial role in control- 
ling translational accuracy, and its efficiency can be further modulated by posttranslational 
modifications [64,65]. It is possible that Rps2-mediated noise regulation occurs at the protein 
translation level. More experiments are required to address this issue. 

A recent study has shown that gene expression noise can influence drug resistance and the 
selection of mutations even in the cells derived from a multicellular organism [66]. The regula- 
tory functions of Hmtl are generally well conserved in human cells, but they are executed by 
six orthologs [36,67]. Human orthologs also exist for many Hmt1 substrates, which interact 
respectively with different human type I enzymes [37,68-70]. For example, the human Hmtl 
ortholog protein arginine N-methyltransferase 4 (PRMT4) interacts with the Snf2 ortholog 
Brahma/SWI2-related gene 1 (BRG1) to facilitate the ATPase activity of the entire human 
SWI/SNEF complex [70]. It is likely that the interactions between Hmt] and its substrates repre- 
sent an ancient regulatory network. In support of this supposition, we observed increased 
noise in S. pombe and S. cerevisiae hmt1 mutants, suggesting that Hmt1-mediated buffering 
evolved hundreds of millions of years ago, before the ancestors of these two species diverged. 
Hmt1-mediated noise buffering may also help populations of unicellular S. pombe to survive 
in stressful environments. However, with regard to multicellular organisms, methylation-regu- 
lated noise buffering may influence their developmental plasticity, and this interesting topic 
awaits further study. 


Materials and methods 
Yeast strains and genetic procedures 


All S. cerevisiae strains used in this study were derived from the W303 strain (Jeu2-3, 112 trp1- 
1 can1-100 ura3-1 ade2-1 his3-11, 15), and the S. pombe strains were derived from the 972 h™ 
strain. Unless otherwise indicated, gene deletion or insertion was based on homologous 
recombination. Yeast cells were transformed using the lithium acetate method [71] or electro- 
poration under 1,800 V/200 (2/25 uF (BTX Gemini SC’, Fisher Scientific, Pittsburgh, PA, 
United States) [72]. 

All fluorescent fusion proteins were constructed by chromosomal in-frame insertion of the 
fluorescent protein tag at the 3’ end of the coding region of target genes. The GFP tags were 
directly amplified from the genomic DNA of corresponding strains in the yeast GFP collection 
[73]. The BFP tag was amplified from the plasmid pFA6a-link-yomTagBFP2-Kan [74]. Prior 
to constructing ADK1-GFP and TDH2 promoter-driven GFP in ancestral and evolved 
TDH2-GFP-carrying lines, the GFP of TDH2-GFP in these strains was removed and replaced 
with a stop codon and the TDH2 3’ untranscribed region (UTR). For TDH2 promoter-driven 
GFP, 1,000 base pairs (bp) upstream of the TDH2 coding region was fused with GFP, and the 
fused fragment was inserted between positions 204886 and 204887 of chromosome I [75]. For 
YFP tagging in S. pombe, the coding region of yVenus [76] without the start codon was fused 
with hphMX6 by two-fragment PCR. To delete S. cerevisiae genes, KanMX4-containing DNA 
fragments for homologous recombination were amplified from the genomic DNA of corre- 
sponding strains in the yeast deletion collection [77]. To construct a hypomorphic mutant of 
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RPS2, we directly amplified the rps2-DAmpP allele from the yeast DAmP diploid collection [78] 
and used it to replace the native gene. To delete rmt1 in S. pombe, a KanMX6-containing DNA 
fragment with homologous flanking regions was generated and used to replace the whole 
gene. 

All single nucleotide polymorphism (SNP) mutant reconstitution strains were constructed 
using the CRISPR/Cas9 system [79,80]. Briefly, host cells were transformed with the Cas9 
expression plasmid, and then the transformant was introduced with a DNA fragment encod- 
ing gRNA, a linearized vector, and mutation-containing donor DNA fragments. After trans- 
formation, single colonies were streaked out on 1% yeast extract, 2% peptone, and 2% dextrose 
(YPD) plates to purify the CRISPR transformants. Genomic DNA of the resulting single colo- 
nies was isolated and examined initially by allele-specific PCR [81] and then by Sanger 
sequencing. The plasmids for Cas9 and gRNA expression are listed in S3 Table. 


Experimental evolution 


For the evolution experiment, eight genes (ADK1, APA1, PCM1, RPL4B, SAM4, TDH2, TPD3, 
and TYS1) were selected and fused with GFP to generate evolving strains, each carrying an 
individual reporter gene. These reporter genes were chosen based on the following criteria: (1) 
each gene belongs to a distinct cellular pathway, (2) the expression level of GFP fusion proteins 
is constant during the cell cycle and is sufficiently high to be detectable by flow cytometry, and 
(3) the reporter proteins are evenly localized in the nucleus or cytosol to avoid misinterpreta- 
tion due to dynamic organelles. 

Because of the complexity of our selection regime, only one evolving line for each reporter 
gene was established. In each selection cycle, 1 x 10° cells were treated with 2.8% EMS (see the 
EMS mutagenesis section below for details) to increase the genetic diversity of the cell popula- 
tion. This mutagenic treatment is crucial, as we initially ran a pilot experiment using a similar 
selection regime but lacking the mutagen treatment and did not observe any obvious increase 
in noise after 70 cycles of selection. However, our EMS treatments also significantly increased 
the number of mitochondria-defective cells, which are known to increase population heteroge- 
neity. To constrain the population of mitochondria-defective cells, we grew cells in nonfer- 
mentable 1% yeast extract, 2% peptone, and 2% glycerol (YPG) medium at 28 °C after EMS 
treatments. After 12 h of growth in YPG, cells were sorted to select 5,000 cells from the top (or 
bottom) 5% of the total population in terms of GFP intensity (see the Fluorescence-activated 
cell sorting section below for details). These cells were grown in 3 ml YPG for approximately 
36 h to reach ODgo0 = 1 before proceeding to the next cycle of selection. We used alternating 
selection between the top 5% and bottom 5% of total populations throughout the evolution 
experiment. The effective population size was estimated to be 1.33 x 10° cells using the formula 
2/N. = 1/(Noi X g) + 1/(Noz X g); in which Np is the initial population size and g is the number 
of generations during each growth period [82]. 

After 35 cycles of selection, five individual clones were isolated from each evolved popula- 
tion, and their GFP noise levels were measured. The clone with the greatest increase in noise 
without exhibiting a decrease in the mean intensity of GFP or increased variation in cell size 
was selected for further genetic analysis. 


EMS mutagenesis 


Mutagenesis was performed according to a previously published protocol by which mutation 
rates can be increased without inducing considerable cell death [83]. Briefly, we washed 

1 x 10° cells with sterile water once and with 100 ul of phosphate buffer (0.1 M Na,HPO, [pH 
7.0]) once and then resuspended them in 90 ul of phosphate buffer. We added 90 ul of EMS- 
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containing phosphate buffer to the cell solution, resulting in a final concentration of 2.8% 
EMS. The solution was maintained under constant shaking at room temperature for 30 min, 
before stopping the reaction by adding 50 ul of 25% sodium thiosulfate. After washing with 
sterile water, the cells were transferred to 10 ml YPG. The survival rate of wild-type cells after 
EMS treatment ranged from 60% to 80%. 


Fluorescence-activated cell sorting and flow cytometry 


Cells were suspended in filtered phosphate-buffered saline (PBS) (137 mM NaCl, 2.7 mM KCl, 
10 mM Na,HPO,, and 1.8 mM KH,PO,) and sampled using a BD FACSJazz machine (BD 
Bioscience, Franklin Lakes, NJ, USA) at an event rate of <3,000 cells/s. GFP readouts were 
acquired with bandpass filters of 513/17 nm using laser excitation at 488 nm. For YFP, we used 
laser excitation at 488 nm with a 542/27-nm filter, and for BFP, we used a 450/50-nm filter and 
laser excitation at 405 nm. 

To eliminate interference from small particles in the fluidic system, we excluded particles 
with forward scatter (FSC) signals < 2,570 (as a trigger threshold). We established three hierar- 
chical gates. First, to avoid cell aggregates, we collected cells that had relatively constant signals 
of trigger pulse width along FSC signals. Second, to eliminate possible cell aggregates, we col- 
lected cells that had relatively constant signals of FSC-W along FSC signals. Third, to avoid cell 
debris and abnormally large cells, we collected cells having FSC and side scatter (SSC) signals 
both ranking within 5%-95% of the population. 

For alternating selection in experimental evolution, we applied the trigger threshold and 
the first and third gates. Based on the distribution of GFP levels, we defined the highest or low- 
est 5% of the total gated population and collected 5,000 cells. The collected cells were then 
grown in 3 ml YPG and used for the next run of selection. For fitness assays, we applied the 
trigger threshold and all three gates. Based on the distribution of GFP levels, we defined the 
high and low subpopulations from the total gated population and collected individual cells. A 
1.0 single-drop sorting mode was employed. The viability of the collected cells was then mea- 
sured under different conditions. 


Noise and signal distribution measurement 


We measured the noise of fluorescent protein signals by the Fano factor (a ratio of variance to 
the mean, %) using more than 5,000 early log-phase cells. This measurement is characterized 
by having less interference from the mean and is more sensitive to variation-driven increased 
noise [12,84]. We used the same trigger threshold and the first two hierarchical gates as 
employed for cell sorting in our noise measurement. For the third gate, we used contour plots 
of FSC and SSC signals to collect cells constituting 60% of the second-gated population to 
ensure more homogenous cell size and cell physiology. Fluorescence readouts of the entire 
population were log-transformed and used to calculate the noise. 


Fitness assays under different growth conditions 


To measure HO; resistance, Tdh2-GFP-carrying cells were grown in YPD for 5 d to ensure 
that the cells had entered stationary phase. More than 100 cells from the top 30% or bottom 
30% of the gated population in terms of their Tdh2-GFP intensity were sorted and spotted 
onto plates with or without 4.4 mM H,03. The survival rates were measured after 5 d of incu- 
bation at 28 °C. 

For heat-resistance assays, cells from the top 10% or bottom 10% in Tdh2-GFP intensity 
were sorted into microcentrifuge tubes. The collected cells were divided into two parts. One 
part was immediately placed in a PCR machine to perform heat ramping from 30 °C to 56 °C 
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for 20 min as a heat stress [85]. The other part was placed in fresh YPD media for 2.5 h and 
then subjected to the same heat treatment. Cell survival rates were measured as colony-form- 
ing units after 3 d of incubation at 28 °C. 

As a rebudding assay, more than 100 unsorted stationary-phase cells were placed on agarose 
pads and monitored by time-lapse microscopy at 15-min intervals for at least 12 h [86]. Only 
unbudded cells at the beginning of the recording period were monitored. 


Western blotting 


Log-phase cells were lysed using NaOH lysis [87]. Proteins in the lysates were separated with 
SDS-PAGE and transferred to a PVDF membrane. Mouse anti-GFP antibody (1:4,000) (#sc- 
9996, Santa Cruz Biotechnology, Dallas, TX, USA) was used to detect Tdh2-GFP. Rabbit anti- 
G6PDH antibody (1:4,000) (#49521, Sigma-Aldrich, St. Louis, MO, USA) was used to detect 
G6PDH, which served as an internal control. Mouse anti-methylated arginine antibody 
(mab0002-P, Covalab, Villeurbanne, France) was used to detect proteins with methylated argi- 
nine. MultiMab rabbit monoclonal mix antibodies against an asymmetric di-methyl arginine 
motif (#13522, Cell Signaling Technology, Danvers, MA, USA) were used to detect proteins 
with asymmetric di-methylated arginine. 


Microscopy 


Microscopy was conducted using a 60x objective lens and an ImageXpress Micro XL system 
(Molecular Device, Sunnyvale, CA, USA). 


Fl segregant analysis 


The evolved line was mated with the ancestral line, and the resultant diploid cells were induced 
to sporulate. Sporulated culture was harvested into a microcentrifuge tube and treated with 0.5 
mg/ml Zymolyase-100T (Nacalai Tesque, Kyoto, Japan) in 1 M sorbitol at 28 °C for 2 h to 
remove the ascal wall. Cells were then treated with 2% SDS at 28 °C for 10 min to kill unsporu- 
lated diploid cells, before washing with sterile water and vortexing vigorously to attach individ- 
ual spores to the tube wall. We then added 0.01% Triton X-100 solution to the tube, before 
vigorous sonication to detach spores from the tube wall and to separate spore clusters. The sus- 
pension was diluted and spread on YPD plates to isolate F1 haploid segregants. 

We conducted a total of three runs of noise measurement to identify “evolved-like” and 
“ancestral-like” F1 segregants. The progenies ranking within the top 20% or bottom 20% of the 
Tdh2-GFP noise level without changing the mean intensity (i.e., within three standard devia- 
tions of the control) were selected for the next run of noise measurement. We started with 360 
F1 progeny of confirmed ploidy in the first round and obtained 16 “evolved-like” segregants 
and 20 “ancestral-like” segregants after finishing the third run. These cells were subjected to 
whole-genome sequencing analysis. 


Whole-genome sequencing analysis 


Yeast cells were resuspended in 200 ul of lysis buffer (2% Triton X-100, 1% SDS, 100 mM 
NaCl, 10 mM Tris-HCl [pH 8.0], and 1 mM Na,EDTA). Glass beads (0.5 mm in diameter) 
that amounted to the volume of the cell pellet were added, followed by the addition of 200 ul of 
phenol:chloroform:isoamyl alcohol (PCIA, 25:24:1 [pH 8.0]). The mixture was vortexed for 5 
min. We then added 200 ul of TE buffer (10 mM Tris-Cl; 1 mM Na,EDTA [pH 8.0]), and the 
mixture was centrifuged for 5 min at maximum speed. Only 350 ul of the solution in the aque- 
ous layer was extracted and subjected to EtOH precipitation. The pellet was air-dried and then 
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dissolved in 400 ul of TE buffer containing 75 ug/ml RNase A. The mixture was then incubated 
at 37 °C for 5 min. We then added 200 ul of PCIA, and we inverted the tubes to mix. EtOH 
precipitation was then performed at room temperature under constant shaking for 2 h. The 
pellet was air-dried and then dissolved in 100 tl of sterile double-distilled water. The ratio of 
OD 609 to OD 29 of the purified DNA ranged from 1.8 to 2.1. Concentrations of gnomic DNA 
were determined using a Qubit dsDNA BR assay kit (ThermoFisher Scientific). 

Equal amounts of genomic DNA from individual “evolved-like” or “ancestral-like” F1 seg- 
regants were pooled. These two DNA pools were sequenced using an Illumina NextSeq500 sys- 
tem (Illumina, San Diego, CA, USA) with 150-bp paired-end reads from 350-bp libraries. At 
least 57x coverage was achieved for sequencing ancestral and evolved clones, and we obtained 
150-200x coverage for the pooled segregants. Sequence results were analyzed using CLC 
Genomics Workbench 9.1 with default settings for read import, read trimming, read align- 
ment, duplicate-read removal, local realignment, SNP calling (with the slight modification of 
using a 10% frequency cut-off in the “ancestral-like” pool to retrieve as many mutations as pos- 
sible, since the default is 35%), and SNP annotation. We used the reference genome of S288C 
(version R64-2-1) from The Saccharomyces Genome Database for SNP calling because it is the 
best-annotated. 


Identifying candidate causal SNPs responsible for increased noise 


By using the CLC Genomics Workbench 9.1 function “filter by control read” with a criterion 
of fewer than 2 control reads, we could eliminate SNPs existing in both ancestral and evolved 
lines (S1 Table). The resulting list of evolved SNPs was then cross-referenced with the list of 
SNPs from the “evolved-like” pool. We anticipated that candidate causal SNPs responsible for 
increased noise would be enriched in the “evolved-like” but not in the “ancestral-like” pools. 
By simulating sequencing results (using Python 2.7) from 150x coverage, we determined the 
frequencies of SNPs that were randomly segregated into two pools of 10 or 20 individuals to 
represent SNP enrichment thresholds for increased noise [88]. We identified SNPs with a 
frequency > 70% in the “evolved-like” pool and with a frequency < 32% in the “ancestral-like” 
pool at a 90% confidence interval as being correlated with increased noise (S2 Table). The 
thresholds are represented by 90% confidence intervals for simulations with 10 individuals 
and at 95% for those with 20 individuals. The script used to determine the threshold of allele 
frequency for the causal mutations was uploaded to GitHub under a repository “imst0715” 
and can be found by searching “Threshold-for-bulk-segregant-analysis.” 


Histone ChIP combined with Q-PCR 


Log-phase cells (approximately 30 OD¢o9) were harvested and cross-linked by 1% formalde- 
hyde at 30 °C for 30 min. After quenching with 125 mM glycine, the cell suspension was lysed 
by three cycles of 5 min beating and 1 min cooling at 4 °C of 0.5-mm glass beads in FA buffer 
(50 mM HEPES-KOH [pH 7.5], 140 mM NaCl, 1% Triton X-100, 0.1% SDS, and 0.1% Deoxy- 
cholate Na salt) supplemented with 1x protease inhibitor cocktail (Protease Inhibitor Cocktail 
Set IV in DMSO; Merck, 539136) and 1 mM PMSF. We made a hole in the bottom of the 
lysate-containing tube to allow the lysate to flow through into a collection tube upon centrifu- 
gation at 500g for 3 min at 4 °C. The chromatin fraction was pelleted down by centrifugation 
at 12,000g for 15 min at 4 °C and washed twice with FA buffer. We transferred all of the chro- 
matin suspension (2 ml) into a 15-ml centrifugation tube (BIOFIL), avoiding bubbles. We 
sheared the chromatin using a Bioruptor (Diagenode, Denville, NJ, USA), with 5 cycles of 15 
min (30 s on and 30 s off per minute) at “High intensity” mode at 4 °C. Ice-cold water was 
resupplied in the Bioruptor tank after each cycle to maintain the low temperature. The sheared 
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chromatin was cleared by centrifugation at 12,000g for 15 min at 4 °C. One-fifteenth of the 
cleared chromatin was used as input and stored at —80 °C. The remaining cleared chromatin 
was transferred into a tube containing preincubated Dynabeads Protein A (#10002D, Invitro- 
gen, Waltham, MA, USA) with histone 3 antibody (ab#1791, Abcam), followed by end-over- 
end mixing overnight at 4 °C. The beads were anchored with DynaMag-2 Magnet (Thermo- 
Fisher Scientific, Waltham, MA, USA) to remove unbound molecules and were then sequen- 
tially washed with 1 ml of FA buffer, 500 mM NaCl-containing FA buffer, DOC buffer (10 
mM Tris-HCl, 1 mM Na,EDTA [pH 8.0], 250 mM LiCl, 0.5% NP-40, 0.5% Deoxycholate Na 
salt), and TE buffer (10 mM Tris-Cl; 1 mM Na,EDTA [pH 8.0]). The bound molecules were 
eluted by adding TES buffer (10 mM Tris-Cl; 1 mM Na,EDTA [pH 8.0], and 1% SDS) and 
incubating at 65 °C for 20 min, followed by a second elution with TE buffer (65 °C for 10 
min). The two eluents were combined. The eluent and thawed input sample were incubated in 
0.125 ug/ml RNase A at 37 °C for 30 min, followed by Proteinase K treatment (2 mg/ml Pro- 
teinase K at 42 °C for 1 h). The samples were de-cross-linked at 65 °C overnight. DNA was 
purified using a QIAquick DNA Purification Kit (Qiagen, Hilden, Germany) according to the 
manufacturer’s instructions, with the modification of a two-cycle wash step using Qiagen PE 
buffer. DNA (0.25 ng) was then subjected to Q-PCR with the primers for nucleosome-occu- 
pied and nucleosome-depleted regions of the promoter and coding regions for TDH2 and 
GLY1 (S3 Table) [89]. The amplicons for TDH2 and GLY] regions were normalized to the 
amplicon of a well-defined nucleosome-occupied region at GAL1 promoter [90]. The nucleo- 
some occupancy was determined by a ratio of the normalized amplicons from the ChIP sam- 
ples to those from the input samples. 


Q-PCR of mRNA under different growth conditions 


For stress conditions, log-phase cells grown in YPD media were divided into four parts. One 
part was continuously grown in YPD as an untreated sample, and the remaining three parts 
were treated with 0.375 M KCI for 20 min or 0.4 mM H,0, for 20 min or subjected to 42 °C 
for 30 min. For different growth states, total RNA was extracted from the same batch of cul- 
tures when cells were in log phase and entering diauxic shift. To investigate the activity of 
SWI/SNF chromatin-remodeling complexes, cells were first grown to log phase in Serine- 
depleted complete synthetic mixture (CSM-Serine). Then, L-Serine was added to a concentra- 
tion of 0.1% to induce CHA1 expression, before harvesting cultures 0, 5, and 15 min after 
induction, with the addition of 0.05% NaN; to stop transcription [52]. 

RNA from 5-10 ODgop of cells was extracted and quantified according to a previous report 
[72] with modifications. Cells were harvested at 4 °C. RNA quality was examined using an Agi- 
lent RNA 6000 Nano kit with an Agilent 2100 Bioanalyzer (Agilent Technologies). RNA (2 ug) 
was reverse-transcribed with 0.5 ug oligo(dT): primer using an Applied Biosystem High- 
Capacity cDNA Reverse Transcription kit (ThermoFisher Scientific). The reaction products 
(1 pl) were then used for Q-PCR in an Applied Biosystem 7500 Fast Real-Time PCR system 
(ThermoFisher Scientific) with gene-specific primers (S3 Table). 


Statistical analyses 


All statistics were performed using R language (http://www.r-project.org/). 


Supporting information 


$1 Fig. Sorted subpopulations of cells exhibiting low or high Tdh2-GFP are probably not 
genetically distinct. Cells were sorted according to their Tdh2-GFP levels into low (red circle) 
and high (blue circle) subpopulations, each of which constituted 10% of the whole population 
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(gray). These cells were grown in YPD for the indicated periods of time and analyzed by flow 
cytometry. The distributions of the low (red) and high (blue) Tdh2-GFP subpopulations are 
superimposed with that of the parental population (gray outline). The vertical dashed lines 
indicate the mean signal intensity of the parental population. Data associated with this figure 
can be found in the supplemental data file (S1 Data). GFP, green fluorescent protein; YPD, 1% 
yeast extract, 2% peptone, and 2% dextrose. 

(EPS) 


$2 Fig. Experimental evolution alters the Tdh2-GFP signal distribution without affecting 
protein integrity or subcellular localization. (A) Increased noise in the evolved Tdh2-GFP- 
carrying line is due to a flattened distribution rather than a bimodal one. Data associated with 
this figure can be found in the supplemental data file (S1 Data). (B) Full-length Tdh2-GFP is 
retained after evolution. Western blots were hybridized using mouse anti-GFP antibody 
(1:4,000) (the upper blot) and rabbit anti-G6PDH antibody (1:4,000) (the lower blot). The 
original image of the blot can be found in the supplemental file (S1 Raw Images). (C) Protein 
localization of Tdh2-GFP is not altered by our experimental evolution approach. Cells were 
imaged using a 60x objective under the FITC channel (the upper panel) or bright field (the 
lower panel). The scale bar is 5 um. FITC, fluorescein isothiocyanate; G6PDH, glucose-6-phos- 
phate dehydrogenase; GFP, green fluorescent protein. 

(EPS) 


S3 Fig. Tdh2-GFP signal in F1 progeny selected for whole-genome sequencing. For bulk 
segregant analysis, 360 Fl progeny were derived by backcrossing the evolved clone to the 
ancestral clone. The noise and mean of Tdh2-GFP signal in individual progeny were analyzed, 
and a total of three runs of noise measurement were conducted to identify “evolved-like” and 
“ancestral-like” F1 segregants. Data from the final run of analysis are shown here. Genomic 
DNA of the “evolved-like” and “ancestral-like” Fl progeny was extracted and respectively 
pooled for whole-genome sequencing. Data associated with this figure can be found in the sup- 
plemental data file (S1 Data). GFP, green fluorescent protein. 

(EPS) 


$4 Fig. The hmt1-G70D mutation phenocopies the loss-of-function mutation. (A) The 
identified G70D mutation of Hmtl is located in a highly conserved methyltransferase motif. 
An alignment of the primary sequences of the conserved motif is shown for various methyl- 
transferases from budding yeast, fission yeast, human, and bacteria. Residues shared with S. 
cerevisiae Hmtl are labeled in yellow, and the mutated glycine residue observed in the evolved 
clone is indicated by an arrowhead. (B) Both hmt1-G70D and Hmtl deletion mutants exhibit a 
similar level of increased Tdh2-GFP noise (one-sided Wilcoxon rank-sum test, n = 5-10; 

p = 0.0013 for hmt1-G70D, p = 0.0013 for hmt1A). Data associated with this figure can be 
found in the supplemental data file (S1 Data). (C) The G70D mutation of Hmtl results in 
defective methylation. The patterns of asymmetric dimethylation (left) and panmethylation 
(right) of arginine in whole-cell lysates of hmt1-G70D mutant cells did not differ from those of 
deletion mutants (hmt1A) but differed from those of wild-type cells. The arrowheads indicate 
the signals differentially detected between methylation-competent and methylation-deficient 
strains. The original image of the blot can be found in the supplemental file (S1 Raw Images). 
GFP, green fluorescent protein. 

(EPS) 


$5 Fig. Tdh2-GFP noise significantly increases shortly after stress treatments. Tdh2-GFP 
noise increases under nonoptimal growth conditions (one-sided Wilcoxon rank-sum test, 
n = 4; p = 0.015 for heat stress, p = 0.015 for oxidative stress). Noise was measured after cells 
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were treated with the indicated stress for 20-30 min. The difference between untreated and 
treated cells became more obvious after 2 h (Fig 6B), which probably reflects the time it takes 
for cells to alter the abundance of Tdh2-GFP protein. The median value of replicates is indi- 
cated with horizontal solid lines among groups of data points. *p < 0.05. Data associated with 
this figure can be found in the supplemental data file (S1 Data). GFP, green fluorescent pro- 
tein. 

(EPS) 


S1 Table. Mutations in the evolved TDH2-GFP-carrying strain. GFP, green fluorescent pro- 
tein. 
(XLSX) 


$2 Table. Identities and frequencies of candidate mutations responsible for increased 
noise, and identified mutations in the ancestral-like and evolved-like F1 progeny pools. 
(XLSX) 


$3 Table. Plasmids and primers used in this study. 
(XLSX) 


S1 Data. Data files related to Figs 1B, 2A-2C, 3A, 3B, 4A-4C, 5A-5E, and 6A-6D, as well 
as S1, S2A, S3, S4B, and S5 Figs. 
(XLSX) 


$1 Raw Images. The original images of the blots shown in S2B and S4C Figs. 
(PDF) 


Acknowledgments 


We thank members of the Leu lab for helpful discussion and comments on the manuscript. 
We thank Shao-Win Wang for the S. pombe strains. We also thank John O’Brien for manu- 
script editing and Po-Hsiang Hung and the IMB Bioinformatics, Genomics, and Imaging 
cores for technical assistance. 


Author Contributions 

Conceptualization: Jun-Yi Leu. 

Formal analysis: Shu-Ting You, Jun-Yi Leu. 
Funding acquisition: Jun-Yi Leu. 

Investigation: Shu-Ting You, Yu-Ting Jhou, Jun-Yi Leu. 
Methodology: Cheng-Fu Kao. 

Project administration: Jun-Yi Leu. 

Resources: Jun-Yi Leu. 

Supervision: Jun-Yi Leu. 

Validation: Shu-Ting You. 

Writing - original draft: Shu-Ting You, Jun-Yi Leu. 


References 


1. Elowitz MB, Levine AJ, Siggia ED, Swain PS. Stochastic gene expression in a single cell. Science. 
2002; 297(5584):1 183-6. https://doi.org/10.1126/science. 1070919 PMID: 12183631. 





PLOS Biology | htips://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 21/26 


©-PLOS | BIOLOGY 


Experimental evolution of noise regulation 





10. 


11. 


12. 


13. 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23. 


24. 


Spudich JL, Koshland DE. Non-Genetic Individuality—Chance in Single Cell. Nature. 1976; 262 
(5568):467—71. https://doi.org/10.1038/262467a0 PMID: 958399 


Suda T, Suda J, Ogawa M. Single-cell origin of mouse hemopoietic colonies expressing multiple line- 
ages in variable combinations. Proc Natl Acad Sci U S A. 1983; 80(21):6689—-93. hitps://doi.org/10. 
1073/pnas.80.21.6689 PMID: 6579554. 


Losick R, Desplan C. Stochasticity and cell fate. Science. 2008; 320(5872):65-8. https://doi.org/10. 
1126/science. 1147888 PMID: 18388284. 


Balazsi G, van Oudenaarden A, Collins JJ. Cellular decision making and biological noise: from microbes 
to mammals. Cell. 2011; 144(6):910—25. https://doi.org/10.1016/j.cell.201 1.01.030 PMID: 21414483. 


Raj A, Rifkin SA, Andersen E, van Oudenaarden A. Variability in gene expression underlies incomplete 
penetrance. Nature. 2010; 463(7283):91 3-8. https://doi.org/10.1038/nature08781 PMID: 20164922. 


Burga A, Casanueva MO, Lehner B. Predicting mutation outcome from early stochastic variation in 
genetic interaction partners. Nature. 2011; 480(7376):250-3. https://doi.org/10.1038/nature10665 
PMID: 22158248. 


Acar M, Mettetal JT, van Oudenaarden A. Stochastic switching as a survival strategy in fluctuating envi- 
ronments. Nature Genetics. 2008; 40(4):471-5. https://doi.org/10.1038/ng.110 PMID: 18362885 


Levy SF, Ziv N, Siegal ML. Bet hedging in yeast by heterogeneous, age-correlated expression of a 
stress protectant. PLoS Biol. 2012; 10(5):e1001325. https://doi.org/10.137 1/journal.pbio. 1001325 
PMID: 22589700. 


de Jong IG, Haccou P, Kuipers OP. Bet hedging or not? A guide to proper classification of microbial sur- 
vival strategies. Bioessays. 201 1; 33(3):215—23. https://doi.org/10.1002/bies.201000127 PMID: 
21254151. 


Kaern M, Elston TC, Blake WJ, Collins JJ. Stochasticity in gene expression: from theories to pheno- 
types. Nat Rev Genet. 2005; 6(6):451-64. https://doi.org/10.1038/nrg1615 PMID: 15883588. 


Raser JM, O’Shea EK. Control of stochasticity in eukaryotic gene expression. Science. 2004; 304 
(5678):181 1-4. https://doi.org/10.1126/science.1098641 PMID: 15166317. 


Bar-Even A, Paulsson J, Maheshri N, Carmi M, O’Shea E, Pilpel Y, et al. Noise in protein expression 
scales with natural protein abundance. Nat Genet. 2006; 38(6):636—43. hitps://doi.org/10.1038/ng1807 
PMID: 16715097. 


Taniguchi Y, Choi PJ, LiGW, Chen H, Babu M, Hearn J, et al. Quantifying E. coli proteome and tran- 
scriptome with single-molecule sensitivity in single cells. Science. 2010; 329(5991):533-8. https://doi. 
org/10.1126/science. 1188308 PMID: 20671 182. 


Newman JR, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, et al. Single-cell proteomic 
analysis of S. cerevisiae reveals the architecture of biological noise. Nature. 2006; 441(7095):840-6. 
https://doi.org/10.1038/nature04785 PMID: 16699522. 


Barroso GV, Puzovic N, Dutheil JY. The Evolution of Gene-Specific Transcriptional Noise Is Driven by 
Selection at the Pathway Level. Genetics. 2018; 208(1):173-89. https://doi.org/10.1534/genetics.117. 
300467 PMID: 29097405 


Charlebois DA. Effect and evolution of gene expression noise on the fitness landscape. Phys Rev E. 
2015; 92(2). https://doi.org/10.1103/PhysRevE.92.022713 PMID: 26382438 


Bahar R, Hartmann CH, Rodriguez KA, Denny AD, Busuttil RA, Dolle ME, et al. Increased cell-to-cell 
variation in gene expression in ageing mouse heart. Nature. 2006; 441(7096):101 1-4. https://doi.org/ 
10.1038/nature04844 PMID: 16791200. 


Charlebois DA, Balazsi G, Kaern M. Coherent feedforward transcriptional regulatory motifs enhance 
drug resistance. Phys Rev E. 2014;89(5). https://doi.org/10.1103/PhysRevE.89.052708 PMID: 
25353830 


Hsieh YY, Hung PH, Leu JY. Hsp90 regulates nongenetic variation in response to environmental stress. 
Mol Cell. 2013; 50(1):82—92. https://doi.org/10.1016/j.molcel.2013.01.026 PMID: 23434373. 


Sanchez A, Choubey S, Kondev J. Regulation of noise in gene expression. Annu Rev Biophys. 2013; 
42:469-91. https://doi.org/10.1146/annurev-biophys-083012-130401 PMID: 23527780. 


Battich N, Stoeger T, Pelkmans L. Control of Transcript Variability in Single Mammalian Cells. Cell. 
2015; 163(7):1596—610. https://doi.org/10.1016/j.cell.2015.11.018 PMID: 26687353. 


Ansel J, Bottin H, Rodriguez-Beltran C, Damon C, Nagarajan M, Fehrmann §S, et al. Cell-to-cell stochas- 
tic variation in gene expression is a complex genetic trait. PLoS Genet. 2008; 4(4):e1000049. https:// 
doi.org/10.1371/journal.pgen.1000049 PMID: 18404214. 


Wu §, LiK, Li Y, Zhao T, LiT, Yang YF, et al. Independent regulation of gene expression level and 
noise by histone modifications. PLoS Comput Biol. 2017; 13(6):e1005585. https://doi.org/10.1371/ 
journal.pcbi.1005585 PMID: 28665997. 





PLOS Biology | htips://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 22/26 


@PLOS | BIOLOGY 


Experimental evolution of noise regulation 





25. 


26. 


27. 


28. 


29. 


30. 


31. 


32. 


33. 


34. 


35. 


36. 


37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


Freed NE, Silander OK, Stecher B, Bohm A, Hardt WD, Ackermann M. A simple screen to identify pro- 
moters conferring high levels of phenotypic noise. PLoS Genet. 2008; 4(12):e1000307. https://doi.org/ 
10.1371/journal.pgen.1000307 PMID: 19096504. 


Carey JN, Mettert EL, Roggiani M, Myers KS, Kiley PJ, Goulian M. Regulated Stochasticity in a Bacte- 
rial Signaling Network Permits Tolerance to a Rapid Environmental Change. Cell. 2018; 173(1):196- 
207 e14. https://doi.org/10.1016/j.cell.2018.02.005 PMID: 29502970. 


Grant CM, Quinn KA, Dawes IW. Differential protein S-thiolation of glyceraldehyde-3-phosphate dehy- 
drogenase isoenzymes influences sensitivity to oxidative stress. Mol Cell Biol. 1999; 19(4):2650-6. 
https://doi.org/10.1128/mcb.19.4.2650 PMID: 10082531. 


Cohen D. Optimizing reproduction in a randomly varying environment. J Theor Biol. 1966; 12(1):119- 
29. https://doi.org/10.1016/0022-5193(66)90188-3 PMID: 6015423. 


Venable DL. Bet hedging in a guild of desert annuals. Ecology. 2007; 88(5): 1086-90. https://doi.org/10. 
1890/06-1495 PMID: 17536393 


Graham JK, Smith ML, Simons AM. Experimental evolution of bet hedging under manipulated environ- 
mental uncertainty in Neurospora crassa. Proc Biol Sci. 2014; 281 (1787). https://doi.org/10.1098/rspb. 
2014.0706 PMID: 24870047. 


Stewart-Ornstein J, Weissman JS, El-Samad H. Cellular noise regulons underlie fluctuations in Saccha- 
romyces cerevisiae. Mol Cell. 2012; 45(4):483-93. https://doi.org/10.1016/j.molcel.2011.11.035 PMID: 
22365828. 


Konrad M. Analysis and in vivo disruption of the gene coding for adenylate kinase (ADK1) in the yeast 
Saccharomyces cerevisiae. J Biol Chem. 1988; 263(36):19468-74. PMID: 2848829. 


McNeil JB, McIntosh EM, Taylor BV, Zhang FR, Tang S, Bognar AL. Cloning and molecular characteri- 
zation of three genes, including two genes encoding serine hydroxymethyltransferases, whose inactiva- 
tion is required to render yeast auxotrophic for glycine. J Biol Chem. 1994; 269(12):9155-65. PMID: 
8132653. 


Ahmad I, Rao DN. Functional analysis of conserved motifs in EcoP 151 DNA methyltransferase. J Mol 
Biol. 1996; 259(2):229—40. https://doi.org/10.1006/jmbi.1996.0315 PMID: 8656425. 


McBride AE, Weiss VH, Kim HK, Hogle JM, Silver PA. Analysis of the yeast arginine methyltransferase 
Hmt1p/Rmt1p and its in vivo function—Cofactor binding and substrate interactions. Journal of Biological 
Chemistry. 2000; 275(5):3128-36. https://doi.org/10.1074/jbc.275.5.3128 PMID: 10652296 


Gary JD, Lin WJ, Yang MC, Herschman HR, Clarke S. The predominant protein-arginine methyltrans- 
ferase from Saccharomyces cerevisiae. J Biol Chem. 1996; 271(21):12585—94. https://doi.org/10.1074/ 
jbc.271.21.12585 PMID: 8647869. 


Li HT, Gong T, Zhou Z, Liu YT, Cao X, He Y, et al. Yeast Hmt1 catalyses asymmetric dimethylation of 
histone H3 arginine 2 in vitro. Biochem J. 2015; 467(3):507—-15. https://doi.org/10.1042/BJ20141437 
PMID: 25715670. 


Jackson CA, Yadav N, Min §, Lid, Milliman EJ, Qu J, et al. Proteomic analysis of interactors for yeast 
protein arginine methyltransferase Hmt1 reveals novel substrate and insights into additional biologi- 
cal roles. Proteomics. 2012; 12(22):3304—14. https://doi.org/10.1002/pmic.201200132 PMID: 
22997150 


Lacoste N, Utley RT, Hunter JM, Poirier GG, Cote J. Disruptor of telomeric silencing-1 is a chromatin- 
specific histone H3 methyltransferase. Journal of Biological Chemistry. 2002; 277(34):30421-4. https:// 
doi.org/10.1074/jbc.C200366200 PMID: 12097318 


Yu MC, Lamming DW, Eskin JA, Sinclair DA, Silver PA. The role of protein arginine methylation in the 
formation of silent chromatin. Genes Dev. 2006; 20(23):3249-54. https://doi.org/10.1101/gad.1495206 
PMID: 17158743. 


Wong CM, Tang HM, Kong KY, Wong GW, Qiu H, Jin DY, et al. Yeast arginine methyltransferase 
Hmt1p regulates transcription elongation and termination by methylating NpI3p. Nucleic Acids Res. 
2010; 38(7):2217-28. https://doi.org/10.1093/nar/gkp1133 PMID: 20053728. 


Brandariz-Nunez A, Zeng F, Lam QN, Jin H. Sbp1 modulates the translation of Pab1 mRNA in a poly 
(A)- and RGG-dependent manner. RNA. 2018; 24(1):43-55. htips://doi.org/10.1261/rna.062547.117 
PMID: 28986506. 


Lipson RS, Webb KJ, Clarke SG. Rmt1 catalyzes zinc-finger independent arginine methylation of ribo- 
somal protein Rps2 in Saccharomyces cerevisiae. Biochem Biophys Res Commun. 2010; 391 
(4):1658-62. https://doi.org/10.1016/j.bbre.2009. 12.112 PMID: 20035717. 


Bachand F, Silver PA. PRMT3 is a ribosomal protein methyltransferase that affects the cellular levels of 
ribosomal subunits. Embo Journal. 2004; 23(13):2641-50. https://doi.org/10.1038/sj.emboj.7600265 
PMID: 15175657 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 23/26 


@-PLOS | BIOLOGY 


Experimental evolution of noise regulation 





45. 


46. 


47. 


48. 


49. 


50. 


51. 


52. 


53. 


54. 


55. 


56. 


57. 


58. 


59. 


60. 


61. 


62. 


63. 


64. 


65. 


Poornima G, Shah S, Vignesh V, Parker R, Rajyaguru PI. Arginine methylation promotes translation 
repression activity of elF4G-binding protein, Scd6. Nucleic Acids Res. 2016; 44(19):9358-68. https:// 
doi.org/10.1093/nar/gkw762 PMID: 27613419. 


Henry MF, Silver PA. A novel methyltransferase (Hmt1p) modifies poly(A)+-RNA-binding proteins. Mol 
Cell Biol. 1996; 16(7):3668—78. https://doi.org/10.1128/mcb.16.7.3668 PMID: 8668183. 


Shen EC, Henry MF, Weiss VH, Valentini SR, Silver PA, Lee MS. Arginine methylation facilitates the 
nuclear export of hnRNP proteins. Genes Dev. 1998; 12(5):679-91. https://doi.org/10.1101/gad.12.5. 
679 PMID: 9499403. 


Green DM, Marfatia KA, Crafton EB, Zhang X, Cheng XD, Corbett AH. Nab2p is required for poly(A) 
RNA export in Saccharomyces cerevisiae and is regulated by arginine methylation via Hmt1p. Journal 
of Biological Chemistry. 2002; 277(10):7752-60. https://doi.org/10.1074/jbc.M1 10053200 PMID: 
11779864 


Chen YC, Milliman EJ, Goulet I, Cote J, Jackson CA, Vollbracht JA, et al. Protein arginine methylation 
facilitates cotranscriptional recruitment of pre-mRNA splicing factors. Mol Cell Biol. 2010; 30(21):5245— 
56. https://doi.org/10.1128/MCB.00359-10 PMID: 20823272. 


Laurent BC, Treitel MA, Carlson M. Functional interdependence of the yeast SNF2, SNF5, and SNF6 
proteins in transcriptional activation. Proc Natl Acad Sci U S A. 1991; 88(7):2687-91. https://doi.org/10. 
1073/pnas.88.7.2687 PMID: 1901413. 


Cairns BR, Kim YJ, Sayre MH, Laurent BC, Kornberg RD. A multisubunit complex containing the SWI1/ 
ADR6, SWI2/SNF2, SWI3, SNF5, and SNF6 gene products isolated from yeast. Proc Natl Acad Sci US 
A. 1994; 91(5):1950-4. hitps://doi.org/10.1073/pnas.91.5.1950 PMID: 8127913. 


Ansari SA, Paul E, Sommer S, Lieleg C, He QY, Daly AZ, et al. Mediator, TATA-binding Protein, and 
RNA Polymerase II Contribute to Low Histone Occupancy at Active Gene Promoters in Yeast. Journal 
of Biological Chemistry. 2014; 289(21):14981-95. https://doi.org/10.1074/jbc.M113.529354 PMID: 
24727477 


Alejandro-Osorio AL, Huebert DJ, Porcaro DT, Sonntag ME, Nillasithanukroh S, Lwill J, et al. The his- 
tone deacetylase Rpd3p is required for transient changes in genomic expression in response to stress. 
Genome Biology. 2009; 10(5). https://doi.org/10.1186/gb-2009-10-5-r57 PMID: 19470158 


Worley J, Luo X, Capaldi AP. Inositol pyrophosphates regulate cell growth and the environmental stress 
response by activating the HDAC Rpd3. Cell Rep. 2013; 3(5):1476—82. htips://doi.org/10.1016/j. 
celrep.2013.03.043 PMID: 23643537. 


Taylor JW, Berbee ML. Dating divergences in the Fungal Tree of Life: review and new analyses. Myco- 
logia. 2006; 98(6):838-49. PMID: 17486961. 


Booth IR. Stress and the single cell: intrapopulation diversity is a mechanism to ensure survival upon 
exposure to stress. Int J Food Microbiol. 2002; 78(1—2):19-30. hitps://doi.org/10.1016/s0168-1605(02) 
00239-8 PMID: 12222634. 


Fraser D, Kaern M. A chance at survival: gene expression noise and phenotypic diversification strate- 
gies. Mol Microbiol. 2009; 71(6):1333—40. https://doi.org/10.1111/j.1365-2958.2009.06605.x PMID: 
19220745. 


Veening JW, Smits WK, Kuipers OP. Bistability, epigenetics, and bet-hedging in bacteria. Annu Rev 
Microbiol. 2008; 62:193-210. https://doi.org/10.1146/annurev.micro.62.081307.163002 PMID: 
18537474. 


Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for 
interaction datasets. Nucleic Acids Res. 2006; 34(Database issue):D535-9. https://doi.org/10.1093/ 
nar/gkj109 PMID: 16381927. 


Gopinath RK, You ST, Chien KY, Swamy KB, Yu JS, Schuyler SC, et al. The Hsp90-dependent prote- 
ome is conserved and enriched for hub proteins with high levels of protein-protein connectivity. Genome 
Biol Evol. 2014; 6(10):2851-65. https://doi.org/10.1093/gbe/evu226 PMID: 25316598. 


Albert R, Jeong H, Barabasi AL. Error and attack tolerance of complex networks. Nature. 2000; 
406(6794):378-82. https://doi.org/10.1038/35019019 PMID: 10935628. 


Bergman A, Siegal ML. Evolutionary capacitance as a general feature of complex gene networks. 
Nature. 2003; 424(6948):549-52. https://doi.org/10.1038/nature01765 PMID: 12891357 


Levy SF, Siegal ML. Network hubs buffer environmental variation in Saccharomyces cerevisiae. PLoS 
Biol. 2008; 6(11):e264. https://doi.org/10.1371/journal.pbio.0060264 PMID: 18986213. 


Synetos D, Frantziou CP, Alksne LE. Mutations in yeast ribosomal proteins S28 and S4 affect the accu- 
racy of translation and alter the sensitivity of the ribosomes to paromomycin. Biochim Biophys Acta. 
1996; 1309(1-2):156—66. https://doi.org/10.1016/s0167-4781 (96)00128-5 PMID: 8950190. 


Rother S, Strasser K. The RNA polymerase II CTD kinase Ctk1 functions in translation elongation. 
Genes Dev. 2007; 21(11):1409-21. https://doi.org/10.1101/gad.428407 PMID: 17545469. 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 24/26 


©-PLOS | BIOLOGY 


Experimental evolution of noise regulation 





66. 


67. 


68. 


69. 


70. 


71. 


72. 


73. 


74. 


75. 


76. 


77. 


78. 


79. 


80. 


81. 


82. 


83. 


84. 


85. 


86. 


87. 


Farquhar KS, Charlebois DA, Szenk M, Cohen J, Nevozhay D, Balazsi G. Role of network-mediated 
stochasticity in mammalian drug resistance. Nature Communications. 2019;10. 


Blanc RS, Richard S. Arginine Methylation: The Coming of Age. Mol Cell. 2017; 65(1):8-24. https://doi. 
org/10.1016/j.molcel.2016.11.003 PMID: 28061334. 


Low JKK, Wilkins MR. Protein arginine methylation in Saccharomyces cerevisiae. Febs Journal. 2012; 
279(24):4423-43. https://doi.org/10.1111/febs. 12039 PMID: 23094907 


Bedford MT, Clarke SG. Protein Arginine Methylation in Mammals: Who, What, and Why. Molecular 
Cell. 2009; 33(1):1—-13. https://doi.org/10.1016/j.molcel.2008.12.013 PMID: 19150423 


Xu W, Cho H, Kadam S, Banayo EM, Anderson S, Yates JR 3rd, et al. A methylation-mediator complex 
in hormone signaling. Genes Dev. 2004; 18(2):144-56. https://doi.org/10.1101/gad.1141704 PMID: 
14729568. 


Ito H, Fukuda Y, Murata K, Kimura A. Transformation of intact yeast cells treated with alkali cations. J 
Bacteriol. 1983; 153(1):163-8. PMID: 6336730. 


Hsu PC, Yang CY, Lan CY. Candida albicans Hap43 Is a Repressor Induced under Low-lron Conditions 
and Is Essential for lron-Responsive Transcriptional Regulation and Virulence. Eukaryotic Cell. 2011; 
10(2):207-25. https://doi.org/10.1128/EC.00158-10 PMID: 21131439 


Huh WK, Falvo JV, Gerke LC, Carroll AS, Howson RW, Weissman JS, et al. Global analysis of protein 
localization in budding yeast. Nature. 2003; 425(6959):686—91 .. https://doi.org/10.1038/nature02026 
PMID: 14562095. 


Subach OM, Cranfill PU, Davidson MW, Verkhusha VV. An enhanced monomeric blue fluorescent pro- 
tein with the high chemical stability of the chromophore. PLoS ONE. 2011; 6(12):e28674. https://doi. 
org/10.1371/journal.pone.0028674 PMID: 22174863. 


Metzger BP, Yuan DC, Gruber JD, Duveau F, Wittkopp Pu. Selection on noise constrains variation ina 
eukaryotic promoter. Nature. 2015; 521(7552):344—7. https://doi.org/10.1038/nature14244 PMID: 
25778704. 


Nagai T, lbata K, Park ES, Kubota M, Mikoshiba K, Miyawaki A. A variant of yellow fluorescent protein 
with fast and efficient maturation for cell-biological applications. Nat Biotechnol. 2002; 20(1):87-90. 
https://doi.org/10.1038/nbt0102-87 PMID: 11753368. 


Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B, et al. Functional characteri- 
zation of the S-cerevisiae genome by gene deletion and parallel analysis. Science. 1999; 285 
(5429):901-6. https://doi.org/10.1126/science.285.5429.901 PMID: 10436161 


Breslow DK, Cameron DM, Collins SR, Schuldiner M, Stewart-Ornstein J, Newman HW, et al. A com- 
prehensive strategy enabling high-resolution functional analysis of the yeast genome. Nat Methods. 
2008; 5(8):711-8. https://doi.org/10.1038/nmeth.1234 PMID: 18622397. 


DiCarlo JE, Norville JE, Mali P, Rios X, Aach J, Church GM. Genome engineering in Saccharomyces 
cerevisiae using CRISPR-Cas systems. Nucleic Acids Research. 2013; 41(7):4336—43. hitps://doi.org/ 
10.1093/nar/gkt135 PMID: 23460208 


Horwitz AA, Walter JM, Schubert MG, Kung SH, Hawkins K, Platt DM, et al. Efficient Multiplexed Inte- 
gration of Synergistic Alleles and Metabolic Pathways in Yeasts via CRISPR-Cas. Cell Syst. 2015; 
1(1):88-96. https://doi.org/10.1016/j.cels.2015.02.001 PMID: 27135688. 


Liu J, Huang S, Sun M, Liu S, Liu Y, Wang W, et al. An improved allele-specific PCR primer design 
method for SNP marker analysis and its application. Plant Methods. 2012; 8(1):34. https://doi.org/10. 
1186/1746-481 1-8-34 PMID: 22920499. 


Lenski RE, Rose MR, Simpson SC, Tadler SC. Long-Term Experimental Evolution in Escherichia-Coli. 
1. Adaptation and Divergence during 2,000 Generations. American Naturalist. 1991; 138(6):1315—41. 
https://doi.org/10.1086/285289 


Winston F. EMS and UV mutagenesis in yeast. Curr Protoc Mol Biol. 2008;Chapter 13:Unit 13 3B. 
https://doi.org/10.1002/0471 142727.mb1303bs82 PMID: 18425760. 


Ozbudak EM, Thattai M, Kurtser |, Grossman AD, van Oudenaarden A. Regulation of noise in the 
expression of a single gene. Nat Genet. 2002; 31(1):69-73. https://doi.org/10.1038/ng869 PMID: 
11967532. 


Teng X, Cheng WC, Qi B, Yu TX, Ramachandran K, Boersma MD, et al. Gene-dependent cell death in 
yeast. Cell Death Dis. 2011; 2:e188. https://doi.org/10.1038/cddis.201 1.72 PMID: 21814286. 


Lee HY, Cheng KY, Chao JC, Leu JY. Differentiated cytoplasmic granule formation in quiescent and 
non-quiescent cells upon chronological aging. Microb Cell. 2016; 3(3):109-19. https://doi.org/10.15698/ 
mic2016.03.484 PMID: 28357341 


Kushnirov VV. Rapid and reliable protein extraction from yeast. Yeast. 2000; 16(9):857-60. https://doi. 
org/10.1002/1097-0061 (20000630) 16:9<857::AID-YEA561>3.0.CO0;2-B PMID: 10861908 





PLOS Biology | htips://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 25/26 


@-PLOS | BIOLOGY 


Experimental evolution of noise regulation 





88. 


89. 


90. 


91. 


Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, et al. QTL-seq: rapid mapping of quan- 
titative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 
2013; 74(1):174-83. https://doi.org/10.1111/tpj.12105 PMID: 23289725. 


Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, et al. A high-resolution atlas of nucleosome 
occupancy in yeast. Nat Genet. 2007; 39(10):1235—44. https://doi.org/10.1038/ng2117 PMID: 
17873876. 


Brickner DG, Cajigas |, Fondufe-Mittendorf Y, Anmed S, Lee PC, Widom J, et al. H2A.z-mediated locali- 
zation of genes at the nuclear periphery confers epigenetic memory of previous transcriptional state. 
PLoS Biol. 2007; 5(4):704—16. https://doi.org/10.1371/journal.pbio.0050081 PMID: 17373856 


Keren L, van Dijk D, Weingarten-Gabbay S, Davidi D, Jona G, Weinberger A, et al. Noise in gene 
expression is coupled to growth rate. Genome Research. 2015; 25(12):1893-902. hitps://doi.org/10. 
1101/gr.191635.115 PMID: 26355006 





PLOS Biology | https://doi.org/10.1371/journal.pbio.3000433 October 15, 2019 26/26 


